Casimir forces, surface fluctuations, and thinning of superfluid films 
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Recent experiments on the wetting of 4 He have shown that the film becomes thinner at the A 
transition, and in the superfluid phase. The difference in thickness above and below the transition 
has been attributed to a Casimir interaction which is a consequence of a broken continuous symmetry 
in the bulk superfluid. However, the observed thinning of the film is larger than can be accounted 
by this Casimir force. We show that surface fluctuations give rise to an additional force, similar in 
form but larger in magnitude, which may explain the observations. 
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Quantum fluctuations of the electromagnetic field be- 
tween two conducting plates in vacuum result in long- 
ranged attractive interactions. This effect which has only 
recently experimentally verified by a series of high pre- 
cision measurements [if, was first predicted by Casimir 
m 19480. Thirty years later, Fisher and de Gennes 
noted that the confinement of thermal fluctuations in 
fluids leads to similar long-ranged forces Q. Quite gen- 
erally, geometric restrictions on a fluctuating field con- 
strain the normal modes of fluctuations and result in 
fluctuation-induced or Casimir forces. These forces are 
controlled by correlations in the fluid; when the correla- 
tions are long-ranged, corresponding to massless fields, 
Casimir forces decay with distance as a simple power 
law|l H H 0, S H G3 E3 The overall strength of a 
Casimir interaction is typically universal. That is, it 
depends on symmetries of the fluctuating field, and on 
boundary conditions, but not on microscopic details. 

An important example of a Casimir force associated 
with thermal fluctuations in a condensed matter system 
is found in 4 He films at and near the superfluid phase 
transition|l2j. The finite thickness, d, of the film con- 
strains the fluctuations of the superfluid order param- 
eter, which then mediate a Casimir force. Experimen- 
tal demonstration of this force was reported recently by 
Garcia and Chan^^] (GC). To produce films of various 
thicknesses, a stack of copper electrodes were suspended 
at different heights above bulk liquid helium. The thick- 
ness of the wetting layer on each electrode as a function of 
temperature was monitored to gauge the strength of in- 
teractions with the substrate. Figure^shows the change 
in the film thickness Ac? = d—do, as a function of reduced 
temperature t = T — T\, near the superfluid transition 
point for the capacitor labeled "Cap 1" in Ref. The 
quantity do is the thickness of the film well above the A- 
point. As shown in the figure, there is a perceptible de- 
crease in the thickness of the film at the transition point 
t — 0, followed by a substantial drop below the transi- 
tion. The thinning of the film for t > quantitatively 
confirms the theoretical predictions of attractive Casimir 
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FIG. 1: The thinning of the Helium film Ad versus tem- 
perature, based on the data reported in Refs. [l2l llH . The 
dotted line illustrates the expected thinning of the film due 
to a Casimir force resulting from bulk Goldstone modes. It 
also coincides with the magnitude of the thinning of the film 
at the transition temperature as determined in Ref. 0. The 
dashed line shows the expected change in thickness due to 
both bulk Goldstone modes and surface fluctuations. 



interactions between parallel surfaces in the presence of 
Dirichclct boundary conditions |8j. At this point, no the- 
oretical explanation has been put forth for the relatively 
substantial dip in the film thickness at temperatures close 
to and below the bulk transition temperature. In fact, 
the theoretical magnitude of the Casimir force at t = 
is 50 times less than required to give rise to the observed 
maximum thinning 

Hi- 
Further below the transition point, the superfluid film 

partially recovers its thickness, but still remains notice- 
ably thinner than in the normal fluid. The superfluid 
state, in which a continuous symmetry is broken, sup- 
ports Goldstone modes that are not present in the nor- 
mal phase. It is thus reasonable to expect that these long 
wavelength modes are responsible for the differences in 
thickness on the two sides of the transition point. How- 
ever, the magnitude of the fluctuation-induced force as- 
sociated with the Goldstone modes in the bulk of the 
film is too small to account for the observed reduction in 



2 



thickness in the superfluid region away from crit- 

icality). The dotted line in Fig. ^ indicates the expected 
thinning of the film due to the Casimir force resulting 
from the bulk Goldstone modes. 

Here, we explore the role of surface fluctuations on the 
thinning of a superfluid film. In particular, we investigate 
the fluctuation-induced force generated by the flow field 
in a superfluid film that accompanies undulations in the 
fluid-vapor interface. While it may initially appear that 
surface effects should be subdominant to those from the 
bulk, the flow fields (and hence kinetic energy) associated 
with surface deformations are actually quite constrained 
by the film thickness, d. The dependence of undulation 
energies on thickness leads to a contribution to the free 
energy that is proportional to 1/d 2 . This leads to a force 
favoring thinner films, with the same form as that arising 
from the bulk Goldstone modes, and an amplitude that 
is almost twice as large. The dashed line in Fig. ^ shows 
the expected change in the thickness of the film due to 
the combined influences of surface and bulk fluctuation- 
induced forces for a film of thickness 423 A. As is clear 
from the figure, the net effect of "bulk" Goldstone modes 
and surface capillary modes is sufficient to explain the 
diminished breadth of films in the superfluid phase. 

We would like to emphasize that a theoretical calcu- 
lation of forces induced by surface fluctuations in liquid 
films is not without precedent [3. Ilfij], However, we are 
unaware of any instance in which such interactions have 
been detected prior to the experiments of Garcia and 
Chan. The key point is that in most cases, fluctuation- 
induced interactions (due to bulk phonons or surface 
modes) appear only as a small correction to the larger 
van der Waals forces. Traversing the superfluid transi- 
tion, however, provides an instance in which the domi- 
nant atomic forces are unchanged, while the long wave- 
length modes due to continuous symmetry breaking can 
be switched on or off. Thus the change in thickness of the 
wetting film provides an ideal probe of the interactions 
generated by the Goldstone modes. 

Deep in the superfluid state, the magnitude of the com- 
plex order parameter is fixed, but its phase <f> can vary 
with relatively small energy cost. Spatial variations of 
4> are accompanied by the flow of the superfluid compo- 
nent. This connection underlies the two-fluid model of 
superfluidity 0, which leads to the melange of acoustic 
and wavelike excitations supported by this system, in- 
cluding first, second, third and fourth sound 0,0]. For 
our purposes, we decompose the variations of (j> into com- 
ponents arising from bulk and surface modes. In terms of 
the wave- vector k = (k x , k y ,k z = mr/d), the bulk modes 
have the form 



•,(k)=A k e ^+^cos(^) 



(1) 



The cosine in Eq. |Q guarantees that d(j>/dz is equal to 
zero at z — and z = d, so that the flow field associ- 



ated with phase fluctuations does not cause a displace- 
ment of either of the two interfaces. The above mode 
is accompanied by a superfluid velocity v s = HVipb/m, 
and a corresponding kinetic energy (p s /2) J d 3 xv 2 — 
(h 2 p s V/4:m 2 )k 2 \A] si \ 2 , where p s is the superfluid density, 
and m is the mass of a helium atom. 

Integrating over all decompositions of the phase cj> ac- 
cording to Eq. JIJ, we obtain the free energy associated 
with this set of bulk Goldstone modes as 
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Here, A is the total area of the film, and T' corresponds 
to contributions that do not generate a non-trivial de- 
pendence on d. Indeed, the sum over n can be performed 
by standard contour integration techniques 20] ; while the 
dominant term is an extensive contribution to the free en- 
ergy, there is an important correction that scales as 1/d 2 . 
The latter is the fluctuation-induced interaction, which 
leads to a Casimir force per unit area, of 
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The thickness of the adsorbed helium film is deter- 
mined by the competition of several factors, notably the 
loss of gravitational potential energy, and attractions to 
the substrate. The former can be calculated simply from 
the height difference h between the adsorbing plate and 
the bulk liquid, while the latter is due to the van der 
Waals interactions with the substrate^^ El- Fluctua- 
tion induced forces, as in Eq. ©, provide an additional 
component. The film thickness d is thus determined by 
the force balance equation 
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(4) 



The first term on the left hand side is the van der Waals 
interaction, with a leading behavior of 7/d 3 with 7 = 
2600°K A 3 for a film of 4 He on Cu. Retardation effects 
due to the finite speed of light are significant for d of 
the order of di/2 = 193 A, necessitating the inclusion of 
the correction term. The second term on the left hand 
side is the Casimir force, which has the same leading 
behavior as the van der Waals term, with a magnitude set 
b v u = 45.81 A 3 /atom and the dimensionless amplitude 

Unlike the van der Waals interaction, the parameter $ 
is expected to change rapidly at the superfluid transition. 
It is zero in the normal phase, while fluctuations of the 
superfluid order parameter lead to an interaction that is 
'attractive' in the sense of favoring a thinner film. On ap- 
proaching the transition from the normal liquid, the am- 
plitude is a scaling function of d (T — T\) u (y is the ex- 
ponent for the divergence of the correlation length) which 
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was obtained in a two-loop Renormalization Group calcu- 
lation by KrechQ. The predicted thinning of the nearly 
superfluid film closely tracks the observations, e.g. as in 
the portion of Fig. ^ for T >T\. Currently, there are no 
calculations that reproduce the large dip in thickness for 
T < T\fl^. Well below the transition, in the superfluid 
phase, 1} is expected to be a universal constant, such as in 
Eq. JjJJ ■ The maximum amplitude calculated by Krech || 
is coincidentally rather close to the value given in Eq. @ 
coming from the Goldstone modes in the bulk of the film. 
(Monte Carlo simulations with periodic boundary con- 
ditions in three dimensions yield a value of the critical 
amplitude which is roughly twice larger j2||. ) Thus, if 
these modes were the only cause for the thinning of the 
film in the superfluid phase, the height of the film would 
be roughly the same at the critical point and well below 
the A-point. As indicated in Fig.Q] this is not the case. 
More precisely, if h = 0.228 cm (mgh = 10.76 fi°K) the 
equilibrium film thickness can be deduced from Eq. 
to be d — 423 A above the transition temperature. This 
is the thickness of the film for "Cap. 1" in the experi- 
ment of GC. In the case of the bulk Goldstone modes, 
{) = — C(3)/87r, which is too small to explain the reduc- 
tion in thickness below the transition. A straightforward 
calculation indicates that the thinning of the film is con- 
sistent with a force which is at least three times as large 
as the one at the bulk critical point (T — T\). 

To resolve the discrepancy with the experimental ob- 
servation, we now resort to the effect of surface fluctua- 
tions in the superfluid j^. According to Dzyaloshinkskii, 
Lifshitz, and Pitaevskii (DLP), surface fluctuations in flu- 
ids also act as a source for Casimir forces |4|. However, in 
the case of thin liquid films, viscous damping effectively 
clamps the fluid, and there is no indication that such 
forces play any role in the thinning of the helium films. 

But when the film is in the superfluid state, there is no 
viscosity opposing the flow fields that accompany surface 
deformations. To quantify the effect of surface fluctua- 
tions, we consider a film of thickness 



d(R,t)=d + h ? (t) e^, 



(5) 



where q = xq x + yq y is the two-dimensional wave vec- 
tor of the surface distortion, and R is a two-dimensional 
position vector. This deformation is accompanied by a 
distortion in the phase of the superfluid order parameter 
of 



n qsmh(qd) 



(6) 



The form of 4> s is chosen such that the vertical velocity 
v z = (h/m)d z 4> s , is zero at the substrate (z = 0), and 
coincides with the motion of the liquid-vapor interface at 
z = d. Variations of 4> s along the z direction are expo- 



nential in qz, with q = ±./ g| + q 2 . This choice ensures 
that V 2 (j) s = 0, such that the kinetic energy is minimal, 



and that there are no couplings to the bulk modes con- 
sidered earlier. (We note that the velocity potential in 
Eq. reproduces the flow field associated with the third 
sound mode 

The kinetic energy in the flow set up by 4> s is 
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(7) 



Surface deformations are also accompanied by a potential 
energy, and surface tension cost. Including these contri- 
butions, and summing over all surface modes, leads to a 
Hamiltonian 
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In the above equation, we have introduced the conjugate 
momentum 



r - P> h, ; 
q q tanh qd q ' 



(9) 



in the potential energy part p is the mass density, a is 
the surface tension, and / is the net force on the surface 
particles. 

The classical partition function!^] associated with the 
surface modes is easily obtained by integrating over itq 
and hq-, leading to a contribution to the free energy of 
the film of 
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Note that the contribution from the potential energy 
(the integral over hq*) is present both above and below 
the transition, while the kinetic energy term (from in- 
tegration over TTq-) exists only in the superfluid phase. 
In fact, it is only the latter that explicitly depends on 
the thickness of the film through the factor of t&nh(qd): 
the d dependence of the free energy comes entirely from 
(ksT/2) ^,fl n [tanh(gd)], and is independent of various 
material dependent parameters appearing in Eq. (|10|) 28] . 
Taking a derivative with respect to d, we obtain the 
Casimir force resulting from the superfluid flow field in- 
duced by surface fluctuations, as 
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which is nearly twice as large as the "bulk" Goldstone 
mode force in Eq. ©. 

The total Casimir force per unit area due to the surface 
fluctuations and Goldstone modes is 
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Using Eqs. (gj and ifT^I . we find Ad « 2.2A at tempera- 
tures slightly below T\, for a film whose thickness above 
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the transition is d = 423A. The dashed line in Fig. 
illustrates this expected change in thickness of the film, 
which is very close to that observed in the experiments 
of GC. It is important to keep in mind that the interpre- 
tation of wetting experiments on He is complicated by 
the presence of microscopic scratches, dust particles and 
surface roughness. All these have a noticeable effect on 
the estimation of the thickness of the film 0, 0] . 

In conclusion, experiments on helium films adsorbed 
on copper surfaces provide a powerful tool for probing 
fluctuation-induced forces, as the effects due to the su- 
perfluid order parameter are turned on at the transition 
temperature, while the atomic forces are essentially un- 
changed. Fluctuations in the phase of the order param- 
eter are an evident candidate, but can only partially ac- 
count for the observed thinning of the superfluid film. We 
suggest that surface fluctuations provide additional inter- 
actions. In particular, in the superfluid phase the motion 
of the liquid-vapor interface sets up a superfluid velocity 
field that extends throughout the film. The correspond- 
ing fluctuation-induced force has exactly the same (mate- 
rial independent) form as that from the bulk phase fluc- 
tuations, but with an amplitude that is nearly twice as 
large. Considering the statistical and systematic errors 
in the experiments, it is reasonable to conclude that the 
combination of the bulk and surface fluctuation-induced 
forces is sufficient to explain the excess thinning of the 
film in the superfluid region. The relati vely lar ge e rror 
bars on the helium wetting experiments |l2t |l4j, |29j , do 
not mitigate our results, which relate to the difference 
in film thickness above and below the A transition. One 
important extension of this work will be to find the ef- 
fect of surface fluctuations on the Casimir force at — and 
especially immediately below — the superfluid transition. 
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